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ABSTRACT 

We measure the stellar specific angular momentum js = Jsj'^s in four nearby (z ~ 0.1) disk galaxies 
that have stellar masses Mg near the break M* of the galaxy mass function, but look like typical star¬ 
forming disks at 2 : ~ 2 in terms of their low stability (Q ~ 1 ), dumpiness, high ionized gas dispersion 
(40—50 kms“^), high molecular gas fraction (20—30%) and rapid star formation (^ 20 M 0 yr“^). 
Combining high-resolution (Keck-OSIRIS) and large-radius (Gemini-GMOS) spectroscopic maps, only 
available at low 2 :, we discover that these targets have ^3 times less stellar angular momentum than 
typical local spiral galaxies of equal stellar mass and bulge fraction. Theoretical considerations show 
that this deficiency in angular momentum is the main cause of their low stability, while the high gas 
fraction plays a complementary role. Interestingly, the low jg values of our targets are similar to those 
expected in the M*-population at higher 2 ; from the approximate theoretical scaling jg oc (1 -h 
at fixed Mg. This suggests that a change in angular momentum, driven by cosmic expansion, is the 
main cause for the remarkable difference between clumpy M*-disks at high ^ (which likely evolve into 
early-type galaxies) and mass-matched local spirals. 


1 . INTRODUGTION AND MOTIVATION 

In studying the history of the universe, the galaxy pop¬ 
ulation that contributes most of the star formation den¬ 
sity at a given epoch is of central interest. Individual 
galaxies can evolve into and out of this population, but 
remarkably the dominant star-forming population itself 
always remains composed of disks with similar stellar 
masses of Mg = 10^^ —IO^^Mq, at redshifts 2 : < 3 (Karim 
This mass range coincides with the 'knee' 


et al. 2011) 

in the mass function of star-form ing galaxies, c entered at 
M: = lO^^-^-lO^i Mq at ^ < 3 dMuzzin et aL]|2QI3). 


Although the dominant star-forming population always 
consists of M*-galaxies of similar mass, the nature of 
these galaxies changes significantly with z. In modern 
times {z ~ 0 ), most of these galaxies are stable spiral 
disks like the Milky Way. They have low velocity dis¬ 
persion and balance gravity almost entirely by rotation. 
Stars form at modest rates (^ I M 0 yr“^) in thousa nds of 
clouds wi th individual masses around 10^ — 10^ M 0 (Mur- 
ray||20II|). In contrast, galaxies of the same mass in past 
epochs [z ^ 1—3) form stars 10-100 times more rapidly. 
Most of this fast star formation occurs in giant g as clumps 
of 10 ^ —IO^M 0 (Elmegreen & Elmegreen||2005 ), reminis¬ 
cent of merging proto-g alaxies, but now interp reted as 


large Jeans instabilities (Bournand et al. 20071). These 


disks are rich in turbulent gas ( [Genzel et ' aT^ 2bQ8|) and 


yield marginal rotatio nal support and stability (Forster 
Schreiber et al.||2006|). Because of mass growth, these 
clumpy MJ-disks at 2 ; ^ 2 are likely to evol ve into early- 
type galaxies (ETGs, |Tacchella et al.||20I5|). 

The stark morpho-dynamical difterence between early 
and modern star-forming galaxies, in spite of identi¬ 


cal masses, raises the question as to how these galax¬ 
ies ‘know’ what cosmic time they are at. Is their dif¬ 
ference explained by a change in the mass composition, 
such as the establishe d decline in the cold gas fraction 
(Ta cconi et al. _ 2010|) and the mo lecular-to-atomic ra- 
tio (jObreschkow'^fe' Rawlings| 2009), or are these paral- 
lelea by a more fundamental driver? Angular momen¬ 
tum is a promising contender because of its predicted 
2 ;-dependence at fixed galaxy mass: the stellar specific 
angular momentum jg = JgjMg is expected to scale as 
jg (X (1 + 2 :)“^/^ for fixed stellar m ass M g (not individ¬ 
ual galaxies) as discussed in Section 


3.1 This motivates 


the question whether low angular momentum drives or 
contributes to the marginal stability of the clumpy, tur¬ 
bulent M*-disks commonly found at high 2 :. This idea 
is further stimulated by the empirical fact that jg is also 
the dominant driver of morphologic al differences in regu¬ 
lar galaxies of equal mass at 2 : ^ 0 ([Romanowsky fc Eall 


2012 Obreschkow fc Glazebrook||2Ql4| ). 


"Direct angular momentum measurements of clumpy, 
turbulent disks at 2 : > I are, however, challenged by diffi¬ 
culties of conducting deep, high-resolution spectroscopic 
observations at such large cosmic distances. Overcom¬ 
ing this challenge wa s the motivation to assemble the 
‘DYNAMO’ galaxies (Green et al. |20I4), a sample of 
95 galaxies containing some of the most Ha luminous 
(Una > lO^^ergs”^) nearby {z ~ O.I) g alaxies, excludin g 
AGN, in the Sloan Digital Sky Survey (|York et al.|2000|). 
Eollow-up observations (Section]^ revealed that, merg- 
ing objects aside, the Ha luminous DYNAMO objects are 
analogous to typical clumpy disks at 2 : ^ I — 3 (details in 
Section]^. Hence, the DYNAMO sample is an interest- 
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ing laboratory to probe the angular momenta in clumpy 
galaxies, normally too distant for such measurements. 

The objective of this article is to answer two questions: 
(i) what role does angular momentum play in shaping the 
rare clumpy, turbulent disks at low-z, and (ii) what can 
these low-^ studies tell us about the high- 2 : star-forming 
population? These questions will be examined using a pi¬ 
lot sample of four clumpy DYNAMO disks with exquisite 
integral field spectroscopy (ITS) data. Sectionj^describes 
the sample and summarizes the kinematic measurements, 
expanded in the Appendix. Section discusses and in¬ 
terprets the results, synthesized in Section 


2. SAMPLE AND ANALYSIS 

We perform the angular momentum measurements in 
a pilot sample of four DYNAMO galaxies at 2 : 0.1 

(Figure 1 and Table 1), which exhibit Ha luminosities 
in the top 0.1% at their redshifts, but have regular rota¬ 
tion and do not show the typical merger/starburst sig¬ 
natures of Ultra Luminous Infrared Galaxies. The tar¬ 
gets have stellar masses ne ar MN similar to the Milky 
Way {\og{Ms/MQ) = 10.8, |McMillan||2011|), but appear 
to be highly reminiscent 01 typical star-forming galax¬ 
ies at 2 : ^ 2 given the following properties: (i) they are 
composed of giant star-forming c lumps revealed in HST 
Ha-imaging (Bassett et al. 2014), (ii) they exhibit high 
star formation rates (mean of 20 M pyr~^) comparabl e 
to main-sequence galaxies at 2 : > 1 (| Green et al. 2014), 
(hi) they have turbulent (flux weighted gas dispersion 
a = 40—50 kms“^), but rotationally supporte d disks 
[V!(j — 3-5) placed on a Tully-Fisher relation (|Green| 
|et al.), (iv) the disk sizes (i?e ^ 3 kpc) are comparable 
to those of M*-disks at 2 : > 1 (Green et ah), (v) CO(l- 
0) emissio n data uncovers h igh molecular gas fractions 
(20-30%, iFisher et al.|[M4| ). 

Given the strong analogy to high-z turbulent disk 
galaxies, our objects represent an exquisite laboratory for 
studying such clumpy disks at distances 20-times closer 
than possible so far. To our knowledge, the DYNAMO 
sample is unique in that it fulfils all the above properties, 
including ~ Therefore, our sample is a useful 

addition to the well-established class of 4ow-z an alogs 
of high-z objects\ revi ewed by |Glazebrook (2013) and 
Elmegreen et al. (2013). 

Our measurenients of js (details in Appendix A & B) 
rely on high-resolution adaptive optics IFS data from 
Keck-OSIRIS to probe the inner parts and, in two objects 
(G04-1 and G20-2), natural seeing Gemini-GMOS IFS to 
probe the faint outer parts to 2.5i?e- These data give, 
in each pixel, stellar surface densities from the broad¬ 
band infrared/optical continuum and line-of-sight veloci¬ 
ties from Doppler-shifted hydrogen emission lines (Paa 
and H/?). In the case of G04-1 and G20-2, we were 
also able to extract stellar absorption line maps, reveal¬ 
ing that the velocities and dispersions o f the stars are 
consi stent with those of the ionized gas (Bassett et al. 


2014). We can therefore adopt a single rotation curve 
for gas and stars. Where available, data from both spec¬ 
trographs is combined (Appendix A). We note that the 
OSIRIS data confirms the high velocity dispersion found 
by iGreen et al.| (|2014|) (45-50 kms”^ for G04-1, G20- 
2) at a high spatial resolution where the effect of beam 
smearing is negligible (Oliva et al. in prep.). 

Kinematic models are fitted simultaneously to the com¬ 


bined continuum-density and emission line-velocity data, 
making no assumption on the rotation structure other 
than circular rotation in a disk (Appendix B). The fit¬ 
ting routine was successfully verified against degraded 
mock data of 1,000 model-galaxies and the control galaxy 
NGC 3198 (see Appendix G). 

The fits to all targets (Figure 1) are excellent in that 
no global (mono-, di-, or quadrupole) systematics can 
be seen in the density and velocity residuals. Table 1 
lists the fitting parameters and two values for jg, one 
(is,obs) based on the pixels that have IFS data, and one 
{js,tot) that also includes an extrapolation to infinite radii 
assuming an asymptotically flat rotation and exponential 
disk for the outer regions without data. For the two 
galaxies with only inner rotation data, to one effective 
radius Re, from Keck-OSIRIS, the fraction of the total 
specific angular momentum js,tot that is extrapolated is 
about 60%, and the estimated statistical uncertainty of 
js,tot amounts to about 18%. In the two galaxies with 
additional natural seeing Gemini-GMOS data, reaching 
beyond 2 Re, only about 25% of js,tot is extrapolated and 
hence the uncertainty of js,tot reduces to about 11%. This 
shows the utility of such deep IFS data for measuring 
angular momentum. 

In Figure 2, we compare the angular momenta of the 
clumpy targets to typical spiral galaxies at ^ ~ 0 in the 
mass-spin-morphology space, spanned by logM^, logj^, 
and the stellar bulge-to-total mass ratio B/T. The reg¬ 
ular spiral galaxies of The H I Nearby Galaxy Survey 
(THIN GS, [Walter et al.| 2QQ8|) lie near a plane in this 
space (Obreschkow & Glaze brook 2014). In these data, 
the bulge masses in the B/'l' ratios were defined as the 
central exce ss above an exponential den sity profile (Ap¬ 
pendix A of Obreschkow & Glazebrook), irrespective of 
the bulge types (mostly pseudo-bulges). The B/T values 
of our four targets were computed in the same way from 
radial density profiles, extracted from HST r-band (filter 
FR647M) images (for D13-5, G04-1, G20-2) and OSIRIS 
R-band continuum maps (for C22-2). 

For consistency, all stellar masses in Figure 2 were es¬ 
timated from R-band magnitude^ within i nfinite aper¬ 
ture^ usin g a small empirical /^-correction (Glazebrook 
et al!|l9^ ) and a universal R-band mass-to-light raticj^ 


M/Lk — b.5M0/LQ. Using a color-dependent M/Lk 
is more adequate in studies including galaxies w ith sig¬ 
nificant classical b ulges that tend to be redder ( [Fall fc| 


Romanowsky|2013|) . However, in the present comparison 
M, 


a constant M/ suffices: the THINGS and DYNAMO 
galaxies are disk-dominated and t heir continuum maps 


are c onsistent with pseudo-bulges (Kormendy & Fisher 


2008), except for three THINGS objects with small clas¬ 


sical bulges. Overall, the DYNAMO galaxies tend to be 
bluer than typical spirals , but the M/Lk ratio is insen¬ 
sitive (to < 20%) to this ( Bell et ar]|2003 ). 


^ X-magnitudes (Vega zero-point) are from 2MASS for DY¬ 
NAMO and from 3.6p,m Spitzer ma ps with a 2MASS -calibrated 
3.6p,m-to-Ar conversion for THINGS (Leroy et al.|2008|). 

^ For DYNAMO objects, the 2MASS aperture is ~4 scale radii 
containing 99% of the light. For THINGS, the stellar masses in 
(|Leroy et al. |2008|) were calculated using apertures of 1.5r25. By 
extending the aperture to > 5R, we find the fiux to increase by a 
factor 1.2. Th is factor was ap plied to all THINGS galaxies. 


^ Assuing a|Kroupa (|200l|) 
sistency with tJie I HlNUb qj 


, Ini tial Mass Function (IMF) for con- 
ata ( [Leroy et al.|20Q8] >. 
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Table 1 

Properties of the four clumpy target galaxies and the control galaxy NGC 3198. 




Basic properties 

IFS Resolution^ 




Fitted in 

this work^ 



Object 

Redshift 


T ^ 

OSIRIS 

GMOS 

T] 

a 

BIT 

Re 


'dflat 

Js,obs 

js,tot 



log(M0) 

log(ergs“^) 

kpc 

kpc 

deg 

deg 


kpc 

kpc 

km s ^ kpc km s ^ 

kpc km s“ ^ 

G22-2 

0.07II6 

10.4 

41.68 

0.33 

- 

50 

145 

0.12 

2.6 

- 

- 

120 

273 ± 49 

DI3-5 

0.07535 

10.7 

42.10 

0.35 

- 

38 

127 

O.OI 

3.5 

0.5 

202 

356 

828 ± 143 

G04-I 

0.I298I 

II.I 

42.36 

0.57 

1.19 

40 

359 

0.15 

5.1 

0.5 

135 

515 

834± II6 

G20-2 

0.I4II3 

10.7 

42.26 

0.62 

1.28 

46 

187 

0.18 

2.7 

0.6 

II8 

337 

394 ± 32 

NGG 3198 

662 km s“ 

1 10.2 


0.02 

- 

69 

215 

0.00 

6.2 

3.9 

163 

626 

973 ±122 

(reference) 






72^ 

215^ 

0.03-'' 


2.8^ 

150^ 


936 ± 20-'' 


Notes. “ X-band based stellar masses (T^ = 0.5Mq/Lq); ^ Ha luminosity as measured by IPS observations on the AAT/S PIRAL |Green 
eta^20^h ^ FWHM of PSF; ^ Parameters fitted as described in Appendix B; ^ Table 4 of Leroy et al. (2008); ^ Table 1 of Obreschkow & 
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Figure 1. Kinematic fits for the four clumpy targets. Gol 1: emission line fluxes of the Paa (OSIRIS) and (GMOS) lines, as indicated 
in parentheses; is shown in cyan for distinction from Paa. Gols 2&:3: density+velocity maps of data and models, respectively. Lightness 
represents the continuum surface brightness and hue represents the rest-frame line-of-sight velocity of the emission lines. A greyscale is 
used where only continuum data is available. White dashed lines show the best-fit major axes and half-mass ellipses. Gols 4&5: Maps of 
the line-of-sight velocity residuals Av' = v' — "^(nodel stellar density residuals AEg = Eg — Ejnodel,s- 
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3. RESULTS AND DISCUSSION 
3.1. Low js in clumpy galaxies at z ^ 0 

Figure 2 reveals a clear offset of the DYNAMO tar¬ 
gets from the THINGS galaxies, corresponding to a factor 
3.0 ± 0.7 shift in js'. the clumpy targets have three times 
less angular momentum than regular local spiral galaxies 
of comparable Mg and B/T. A likelihood analysis (us¬ 
ing t he Hyper-Fit package of |Robotham fc Obreschkow] 
'2015|), accounting for the statistical uncertainties and 
small number statistics, shows that the offset between 
the two samples is certain at 5.1 standard deviations and 
larger than a factor 2 in js with 97% probability. 

The immediate suggestion of Figure 2 is that the dif¬ 
ference between turbulent, clumpy disks and regular spi¬ 
ral disks is in some sense related to angular momentum. 
The question then becomes one of causality: is the lead¬ 
ing mechanism that (i) low angular momentum causes 
clumps, (ii) clumps cause low angular momentum, or (hi) 
a third quantity causes both? Options (ii) and (iii) are 
hard to substantiate in light of angular momentum be¬ 
ing a conserved quantity that cannot be generated from 
within the galaxy-halo system. Thus, option (iii) would 
mean that external torques somehow cause the giant 
clumps while also reducing the angular momentum - an 
implausible scenario given that our targets show no global 
kinematic asymmetries hinting at such torques. For op¬ 
tion (ii) to work, the star-forming clumps would need to 
drive a mechanism (e.g. via tidal torques or ejected ma¬ 
terial) for transferring angular momentum from the disk 
to the dark halo. However, in m odern high-resolution 
simulations (Genel & et al. 2015) star-forming regions 
preferentially remove low-j gas via super nova-winds, in¬ 
creasing the angular momentum of the disk rather than 
reducing it, in order to make simulated disk galaxies as 
large as observed. In contrast, option (i), i.e. low-j caus¬ 
ing clumps, appears plausible from the following funda¬ 
mental stability considerations. 

The stability of a galactic disk against clump formation 
can be qu antified by the dimensionless Toomre parameter 
( |Tooinre|l964l P 


Q = 


aK 


3.36 GS’ 


( 1 ) 


which depends on the local surface density F, radial ve¬ 
locity dispersion a, and epicyclic frequency k. If Q < 1, 
instabilities can develop with diameters larger than the 
Jeans length Lj and smaller than the maximal stability 
scale I/R set by the differential rotation. If Q > 1, then 
Lj > I/R and no instability can develop: all Jeans modes 
are stabilised by the local centrifugal acceleration. 

In most galactic disks, the surface density and epicyclic 
frequency decrease similarly with radius, suc h that Q 


is constant within 20% (for individual galaxies, Westfall 
|et al.|2014|) across the regions that hold most of the stellar 
mass (at ^ 1 exponential radius) and angular momentum 
(at ^ 2 exponential radii). Within this level of accu¬ 
racy it is therefore possible to introduce a mean Toomre 
parameter Q, defined such that Q > 1 corresponds to 
globally stable disks. The global stability parameter Q 
can be rewritten as Q oc M~^ja upon expressing the 
means of F and n by the total mass M and specific an- 

^ In the case of pure gas disks ctr is commonly expressed via the 
sound-speed Cg, which is well-approximated by the dispersion. 


gular momentum j of the disk (details in Obreschkow 
& Glazebrook 2014). This straightforward substitution 
does not change the underlying physics, but introduces 
a new viewpoint, describing the disk stability in terms 
of fundamental quantities, conserved in isolation. As we 
shall see, this viewpoint allows us to link the disk stability 
to a time-dependent cosmological context. 

So far, Q describes the stability of a single c omponent 
(stellar or gas ) disk. To about 20% accuracy (Romeo & 
Wiegert|[2011[ ) this quantity can be extended to a multn 
component (stars and gas) disk via ( Wang fc Silk||1994| ) 


Q ' = Qd 


Qs 


( 2 ) 


where oc Qg oc Mg ^jgCJg are the partial 

stability terms of the stars and cold gas, resp ectively. 
Relying on the observational finding (see §2 and ~ 


Bassett 


et al. 2014) that stars and gas exhibit comparable velocity 


and dispersion maps, we here adopt the equalities js = 
jg = j and as = (Tg = a. This does not necessarily 
imply that the gas and star s de velop similar instabilities, 
as discussed at the end of ^3.11 


Introducing the gas fraction fg = Mg/{Ms + Afg), the 


global Toomre parameter then reduces to 

Q « - fg)j, a. 


(3) 


Main-sequence spiral galaxies in t he local universe ar e 
stable systems with Q Qiocai = ^ (Westfall et al. |2014| ). 


We now discuss how Q changes relative to Q\q^\ in the 
case of our clumpy turbulent disks. Equation ® shows 
that the high velocity dispersion a of our targets, w hich 


has about twice the amplitude of typical local disks (An- 


dersen et al.||2QQ6 ), would by itself increase the stability 
by a factor 2. 4'he reason is that Lj increases wit h cr, a 


situation that i s also observed in many high-z disks (Gen- 

r 


zel et al. |2QQ8). Thus, the high a can explain the large 


Jeans scale (bisher et al.^ forthcoming), but it cannot si¬ 
multaneously explain why the stability of the disks is low 
enough (i.e. Lr is high enough) for the Jeans modes to 
become manifest. Following Equation ([^, this low sta¬ 
bility despite high a requires high gas fractions or low an¬ 
gular momentum (at fixed Ms). Our targets have both: 
fg = 20—30% (seen in GO emissiorR is a fact or 2-3 


above the typical 10% value of locaUM*-disks (Leroy 
while the js values measured in this work 


et al. 

2008) 

are a 

[actor 


3 lower than in typical local M*-disks. 

Importantly, Equation ^ shows that the factor 1 /3 re¬ 
duction in js has a much larger direct effect on reducing Q 
than the factor 2-3 increase in fg. Thus, low angular mo¬ 
mentum is the dominant driver for the low stability of our 
targets - the main result of this study. Quantitatively, 
the Q-values of our targets are about half those of typical 
local spirals, due to the combination of a factor 0.8 de¬ 
crease in (1 — /g), a factor 0.3 decrease in and a factor 
2 increase in cr, resulting in(5 = 0.8x0.3x2x ~ 1- 

How this low stability feeds back onto cr, via instability- 
driven turbulence and stellar feedback, is a topic of cur¬ 
rent research beyond this study. 

^ The 20—30% refer to molecular material, observed in CO emis¬ 
sion (|Fisher et al.||2014|). An additional atomic (Hi) component 
could increase tliese percentages, but is expected to be situated at 
larger radii, negli gible for the stability of the ste llar component, in 
such dense disks ( [Obreschkow fc Rawlings|2009| >. 
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It is worth discussing some potential caveats of this sta¬ 
bility discussion. Firstly, Equation Q relies on a purely 
mechanical model of a flat d isk. In the case of thick disks, 
the constant 3.36 decreases (|Peng|2001|), while energy dis¬ 
sipation within clumps can increase this constant. In the 
present context, where relative changes in Q matter more 
than its absolute value, these effects only appear at sec¬ 
ond order, unlikely to affect the general results. Secondly, 
we have omitted the fact that the stellar component is 
often less clumpy than the gas - a generic observation 
in most turbulent disks. So the question arises if this 
difference in stellar and gas dumpiness can be explained 
in the context of a single stability parameter Q. While 
this question is beyond the scope of this paper, its an¬ 
swer probably requires a more detailed analysis of the 
dynamics of instabilities. For instance, the evolution of 
a Jeans instability depends on the adiabatic index 7 of 
the fluid: for molecular gas (7 < 4/3), the Jeans mass 
decreases with increasing density and hence any proto¬ 
instability continues its collapse while fragmenting into 
smaller clumps. In turn, the stellar component behaves 
roughly like a monoatomic gas (7 = 5/3), in which the 
Jeans mass increases with increasing density. Hence the 
collapse is stopped making the stellar disk less clumpy. 


3.2. Connection to the clumpy population at z ^ 1 — 3 

Let us now turn to our second objective, the exten¬ 
sion of the discussion to early cosmic times. While our 
four targets are exotic at ^ ~ 0 , their high dispersion, 
dumpiness, star formation and gas fraction are represen¬ 
tative of the main sequence disk population of the same 
stellar mass at z ^ 1 — 3 (Section]^, now thought to 
have formed more star s than merger-driven starbursts 
( Rodighiero et al.|2Qll ). It is clear that the massive star 
formation activity ( 10—100 M 0 yr“^) of this high- 2 : M*- 
population must be fuelled by copious inflows of fresh gas 
(Dekel et al. 2009[ ), which in turn feeds from the inter- 
galactic medium that is ^ 3^ = 27 times denser at 2 : ^ 2 
than today. However, we still require an explanation for 
the low disk stability enabling thi s fresh gas to turn ef - 
ficiently into molecules and stars (Tacconi et al. 2010). 
Is the gas itself the destabilizing agent, or do these high- 
2 : M*-disks exhibit low angular momentum, just as we 
discovered in our four local high- 2 : analogs? 

Measuring the angular momentum in disks at z ^ 2 
is extremely difficult, as their inner rotation structure 
(< 0 ? 1 ) is hard to resolve and their outer parts, holding 
most of the angular momentum, have low surface bright¬ 
ness (owing to the ^ (1 + 2 :)“^ cosmological dimming). 
However, a generic theoretical expectation in ACDM cos¬ 
mology is that the mean relation between mass and an¬ 
gular momentum changes with time, due to cosmic ex¬ 
pansion. To first order, this can be seen in the modei of 
a sphericai haio (viriai mass Mh, r adius Rh, veiocit y W, 
anguiar momentum Jh). Foiiowing |Mo et ah (1998), the 
specific anguiar momentum jh = Jh/Afh = v2XRV can 
be expressed as 

jh = (4) 

where G is the gravitationai c onstant. The spin pa- 
rameter A, here defined as in IBuilock et al.||2QQli is 
appro ximated independent of 2 : (Hetznecker & Burkert 
2006). The Hubbie parameter varies as H[z) = HoR[z), 


where Hq is the iocai Hubbie constant and E{z) = 
(f 2 m(l + -j- in a fiat universe. The function 

Ac{z) is the over-density factor of the viriaiized haio rei- 
ative to the me an density of the univer se, approximated 
anaiyticaiiy as ( Bryan &: Norman||1998 ) 


Ac(^) = 187r^ - 32VL^E-\z) - 390lF-^(z). (5) 

Substituting these anadticai expressions for H{z) and 
Ac( 2 :) into Equation (|^ yieids jh oc M^^^(l + 2 :)“^/^, 
where the 2 :-dependance is an approximation. Assuming 
a z-indepen dent M^/Mh (as verified to 2 : < 4,|Behroozi 
et ai. 2013) and retention of the bar yon specific angu¬ 
lar m omentum during disk formation (Fall & Efstathiou 


1980), a similar ^-scaling applies to js, 

j, OcM2/3(1 + ^)-1/2. (6) 

Modern hydrodynamic galaxy simulations in a cosmo¬ 
logical cont ext, the Magn eticum Pathfinder (Hirsch mann 
et ai.||2014|) and liiustris (jVogeisberger et al.||2014|) sim¬ 
ulations, produce robust angular mome ntum statistics 
dGenel fc et al.||2015| [miu fc et al.| 20 T^ . Both approx- 
imately contLrrn the ^-scaling of Equation (^ for star¬ 
forming M*-disks (priv. comm. S. Genel anoW. Teklu). 

In summary, there is strong theoretical indication that 
the population of M*-disks at z 'zz 2 has lower j^ (by 
a factor ^ 2 ) than the mass-matched disk population at 
z 'zz {). The expected js values of the high- 2 : disks are con¬ 
sistent with to those measured in our local analogs within 
the uncertainties (at 1.3 standard dev iatio ns). Follow¬ 
ing the stability discussion of Section |3.1[ low angular 
momentum must then play a major role in making the 
high- 2 : star-forming disks semi-stable. It thus seems that 
the cosmic growth of angular momentum is a dominant 
driver behind the evolution of the population of star¬ 
forming M/-galaxies from semi-stable, turbulent, clumpy 
to stable, flat, regular spiral disks. We emphasize that 
this conclusion applies to the evolution of the population 
of M*-galaxies, not individual galaxies. A M*-galaxy 
at z ^ 2 does not normally evolve into a M*-galaxy at 
2 : ^ 0, but i nto a much more massive object (van Dokkum 
et al.||2Q13 ) such as an ETG; nor are individual disks ex- 
pected to follow js oc (1 -h z)~^^‘^. 

A potential caveat is that typical z ^ 2 disks of stellar 
mass near M* may have even higher gas fractions (^50%, 
Carilli & Walter 2013) than our objects, which would 
contribute an extra stability factor of 40%. However, 
this is still sub-dominant to the change in j, and may 
arise from it. 


4. GONGLUSIONS 

This article reports the detection of exceptionally low 
angular momentum in four rare examples of clumpy, tur¬ 
bulent star-forming galaxies at 2 : 0.1. The low angu¬ 

lar momenta offer a powerful new explanation for their 
marginal stability, enabling large clumps to form and 
pristine gas to convert rapidly into stars. The key to 
obtaining robust measurements of specific angular mo¬ 
mentum is the combination of IFS data with high spatial 
resolution in the galaxy centers (from Keck-OSIRIS) and 
high surface brightness sensitivity out to large radii (from 
Gemini-GMOS). 

The circumstances through which clumpy low-j disks 
with high gas fractions can exist in the local Universe 
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Clumpy disks: 

• C22-2 

• D13-5 

• G04-1 

• G20-2 

Regular disks: 

• THINGS 


Figure 2. [Interactive 3D-figure in Adobe PDF Viewer. Feel free to use the standalone and animated versions in our Dropbox.] Mass- 
spin-morphology s pace. Dark points represent the 16 typical local spiral galaxies ot the i hlliNGb sample (|Ubreschkow V Giazebrook|2014| 
[Walter et al.|2QQ^ ; colored points are the four clumpy disks of this study - they appear to be deficient in angular momentum jg. Fllipsoids 
represent icr-statistical uncertainties. The jg-values of the DYNAMO objects are the best estimates js,tot (see Appendix B) listed in Table 1. 
The grey plane is best fits to the data accounting for the 3D Gaussian uncertainties. 


remain to be explored, but whatever the origin of these 
objects, it is likely that the physical mechanisms driving 
clump formation and turbulence are the same as at 2 ; ^ 2 
and independent of the galaxies’ origin. By exploiting 
the analogy between our targets and typical star-forming 
galaxies at ^ ~ 2, we find evidence that low angular 
momentum is also a dominant cause of the low stability 
of such high- 2 : galaxies. Hence, a change in angular 
momentum, driven by cosmic expansion, appears to be 
the primary cause for the remarkable difference between 
clumpy M*-disks at high z (which likely evolve into 
ETGs) and mass-matched local spirals. 


We dedicate this publication to the memory Prof Peter McGre¬ 
gor, who passed away during this work. We are indebted to him 
for his outstanding contributions to astrophysics, his leadership of 
Integral Field Spectroscopy in Australia and his inspiration to the 
DYNAMO project. 
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APPENDIX 

A. IFS data and pre-processing 

We observed all four targets with the Keck-OSIRIS integral field spectrograph in Adaptive Optics (AO) mode to 
obtain spectra around the Paa (1875 nm rest-frame) line in 572 x 677 fields. The spatial resolution is 0724 full-width- 
half-maximum (EWHM) of the point spread function (PSE), corresponding to the sub-kpc resolutions listed in Table 1. 
Eor the two more distant galaxies (G04-1 and G20-2) we also obtained deeper lES data around the H/d (486.1 nm rest- 
frame) line using the Gemini-GMOS at a natural seeing of 075 in three fields (one for G04-1, two for G20-2, visible as 
45° rotated rectangles in Eigure 1). 
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The 3D-data cubes (two spatial dimensions, one wavelength dimension) from OSIRIS and CMOS were converted into 
rest-frame line-of-sight velocity maps u'(x) (from Paa, and background-subtracted continuum flux maps Ss(x), on 
a discrete cartesian grid x^ = {xi^yi) with O'/l x O'/l pixels. We assume the stellar surface density to be proportional to 
i;s(x) (hence the subscript). We combine the OSIRIS and CMOS data (for G04-1 and G20-2) by registering the maps 
at 071-pixel accuracy, which suffices given the 075 PSF of GMOS. The GMOS continuum values are scaled to match 
the OSIRIS continuum at the edge of the OSIRIS maps. Merged continuum and velocity maps are then created by 
taking OSIRIS data where available and GMOS data otherwise. In this way, the central parts, where spatial resolution 
is critical, always come from OSIRIS. 


B. Kinematic fits and angular momentum measurement 

The surface density maps line-of-sight velocity maps u'(x) on a discrete grid of 071-pixels x^ = (xi^yi) are 

fitted with a disk model using a customised algorithm: 


1. Fit galaxy centre xq = (^o^^o) by minimizing the convolution of the density-weighted velocity map w = 
and its mirror image, Jd?x u;(x — xo)u;(xo — x). This method is exact for axially symmetric objects. 


2. Fit inclination y and position angle a, which set the projection factor /^^q,(x) = {cos‘^ip + sin^ cpcos 

ip = arg(x — Xq) — a, such that the radius in the galaxy plane is r = /^^q,(x)|x —xo|, and the factor 5f^^Q,(x) = 
/?7,a(x)(cos(/:?sin7^)“^, such that the circular velocity is v{r) = 5f^^cK(x)u'(x). The parameters y and a are fitted 
simultaneously without assuming a rotation or density structure other than a flat, axial disk. This fit is performed 
by minimizing the RMS of the uncertainty-weighted and v residuals in circular rings in the galaxy plane. We 
do not a account for beam-smearing in the velocity field, since the central steep parts of the rotation curve are 
sampled with adaptive optics in all four targets. Explicitly, the maximal relative velocity change per PSF FWHM 
Ax, calculated as max{Vu'(x)}Ax is much smaller than unity in all cases. 


3. Fit a four-parameter disk+bulge density profile Emodei,s(^) = Sdexp(—r/rd) + Ebexp(—r/rb) to Eg using a 
maximum likelihood estimation. Adopting alternative bulge profiles, e.g. with variable Sersic index, makes little 
difference to the quality of these fits and has negligible impact on In the fitting procedure, the model Emodei,s 
is blurred by the PSF, which depends on x due to the combination of adaptive optics and natural seeing. 


4. Fit a two-parameter velocity profile (Boissier et al. 2003) Umodei(^) 
maximum likelihood estimation. 


i'flat(l - exp[-r/rflat]) to v{r) using a 


Using only the first two fitting steps, we compute the directly observed stellar specific angular momentum as 


Js,ohs — 






( 1 ) 


where i goes over all pixels with velocity and continuum data, is the radius in the disk plane and Vi is the circular 
velocity. So far, we only assumed the stellar mass density to be proportional to E^; the proportionality factor turns 
out to be irrelevant for js,ohs as it cancels out in Equation Q. For the same reason, there is no need to correct Eg 
by an inclination factor cos“^ y. The values js,ohs are not converged in the sense that they would increase if data at 
larger radii were available. To obtain the best estimates of the full js, called js,tot^ we extend the sums in Equation Q 
to the domain (to infinite r). Where measurements for E^ or v are unavailable, their values are approximated by 
^modei,s and Vmodeh respectively. 

The best-fit parameters of our targets are listed in Table 1. For G22-2, rflat and v^at are individually almost 
unconstrained, however, their product is constrained at the 15% level, such that the rotation curve is well constrained 
at the radii that most strongly contribute to the angular momentum. We estimate the uncertainties of js,tot by 
propagating the covariance matrix of the fits, taken as the inverse Hessian of the log-likelihood function. For the 
two galaxies with both GMOS and OSIRIS data (G04-1, G20-2), we verified that the extrapolations based solely on 
OSIRIS data {js tot = 958 ± 18%, 397 ± 13%) are consistent with the more accurate extrapolations from combined 
OSIRIS+GMOS’data (j,,tot = 834 ± 14%, 394 ± 8%). 


C. Tests of the fitting algorithm 

We tested the fitting procedure described in Appendix B against 10^ mock datasets of model-galaxies, as illustrated 
in Figure 3. The quality of these mock data is similar to the quality of our targets with OSIRIS+GMOS data. We 
found that our model recovers the true total specific angular momentum js,tot at an accuracy of ^5% with negligible 
systematics. Note that the estimated statistical uncertainties of js,tot in Table 1 are generally larger than 5%, because 
the galaxies are not perfectly axially symmetric disks; hence the fits and extrapolations are somewhat les s certain. 

We also checked the algori thm on the galaxy NGG 3198, a regular M*-disk for which kinematic fits (Leroy et al 


2008) and jg-measurements (Obreschkow & Glazebrook||2014 ) were performed independently at high precision using 
kinematic data reaching out to ^ 10 exponential radii. I^br the present test, the data were heavily truncated and 
smoothed (see Figure 4) to a physical size and resolution similar to the IFS data of the four targets. Applying our 
fitting algorith m to this degraded version of NGG 3198 produces re sults in close agreement with those drawn from the 
reference work ( [Leroy et al.||2008 Obreschkow fc Glazebroo^|2014| ) (bottom of Table 1). 













Mock data ^ M Fitted model 



Velocity residual _ Density residual / ■ 



> 

■A 



Sample index 


Figure 3. LEFT: example of a mock dataset and the corresponding fits with residuals. The mock data represent an inclined model-galaxy 
with exponential surface density profile and differential rotation. The data has been degraded in three ways: (1) it is blurred to a physical 
scale corresponding to the Keck-OSIRIS PSF, (2) signal in some region has been removed (grey and black regions) to mimic sensitivity 
limitations and sensor errors, and (3) Gaussian noise of realistic amplitude has been added. The colorscales are as in Figure 1. RIGHT: 
Relative difference between the recovered value of the total specific angular momentum js,tot and the input value in 10^ random mock 
dataset, similar to the one shown in the left panel. 



Figure 4. Same as Figure 1, but for the heavily degraded local control galaxy NGG 3198 (see Appendix G). 
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